New results on nuclear dependence 
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To understand the reliability of the charmonia suppression as a signature of the Quark- Gluon 
Plasma formation in nucleus-nucleus collisions it is important first to understand the details 
of the production of J/tp and tp' in pA interactions and the difference in the suppression of 
these two states. This report presents the results of the study by the NA50 collaboration of 
the J/ij) and tp' production in pA interactions at 450 GeV beam energy and its dependence on 
rapidity. It is shown that the ip' suffers more suppression than the J/V>, which is consistent 
with a similar observation made at 800 GeV beam energy by the E866/NuSea collaboration. 
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1 Introduction 



The suppression of the production of charmonia states (mainly J/ip and ip') in nucleus — nucleus 
interactions due to colour screening is recognized as one of the most promising QGP signatures^. 
But in real experiments it should be distinguished from the suppression due to the hadronic 
interactions of charmonia with surrounding nuclear matter. This kind of suppression, present 
already in pA interactions, can describe all currently available data except the "anomalous J/ip 
suppression" observed by the NA50 experiment in PbPb collisions @. The deconvolution of 
these two kinds of suppression requires precise quantitative understanding of the interaction of 
different charmonium states with hadronic matter. While a vast amount of data was collected 
for the J/ip, the situation with the ip' is less clear. Published NA50 data showed similar levels 
of suppression for both mesons in pA collisions but a stronger suppression of the ip' in heavier 
colliding systems, like SUB, which was attributed to additional absorption of the ip' due to 
interactions with "comovers" - hadrons produced in the collision. Finally, the E866/NuSea 
collaboration □ observed a stronger suppression for the ip' already in 800 GeV pA interactions. 
The motivation of the present study is to look for similar effects in the NA50 data. 

2 Experimental setup, data samples and analysis 

The main component of the NA50 apparatus is a dimuon spectrometer, covering the 3.0 < yi a f, < 
4.0 rapidity range, which is separated from the target region by a 4.8 m long hadron absorber. 
The dimuon trigger is provided by four scintillator hodoscopes and its efficiency is controlled 
during special runs by two additional hodoscopes. The detailed description of the experiment 
can be found inB. We present here the analysis of data collected in the 1996-2000 period with 
a 450 GeV proton beam interacting with Be, Al, Cu, Ag and W targets. Each data sample was 
obtained with single target, with interaction probability varying from 26 (Al) to 39 (Ag)%. The 
incoming beam intensity was monitored by three argon ionization chambers. For each target 
two qualitatively different data samples were collected: with "low" (LI, 1 — 3 x 10 8 protons per 
2.37 s spill) and "high" (HI, 1 — 3 x 10 9 p/spill) beam intensities. While the LI samples are 
more reliable from the absolute normalization point of view, the HI one provides much better 
statistics, which is necessary for the ip' study. 

The reconstructed dimuons were subjected to —0.5 < y cm < 0.5 and | cos^csl < 0.5 selection 
cuts. For the differential cross section studies, the rapidity range was divided in four equidis- 
tant bins. The mass spectra in each kinematical bin were fitted with curves corresponding to 
individual dimuon sources, obtained from detailed Monte-Carlo simulations. 

Combinatorial background from uncorrelated tt,K — > a decays was estimated from the like- 
sign dimuon samples (using the relation = 2R^J N ++ N , with R being a free parameter 

of the fit, accounting for possible charge correlations in parent hadrons production). The con- 
tribution from the dimuons originated in interactions outside of the target was accounted by 
including into the fit the dimuon spectrum obtained from the special "target-out" runs. 

The systematic errors of the extracted cross sections account for uncertainties in the nor- 
malization: luminosity, trigger and reconstruction efficiencies. 

Since various authors use different parametrizations for the description of charmonia sup- 
pression, we fitted the obtained cross sections by three commonly used models: 

1) The Glauber model, assuming that each charmonium state i is produced in binary nucleon- 
nucleon interactions with cross section <7q and then interacts with surrounding nuclear matter 
with cross section a l abs . This leads to the pA cross section o v a = Vo/v^s / ds* [1 — TA{s)& a bs} A , 
where Ta(s) is the nuclear thickness function at impact vector s. 

2) A simplification of the Glauber model, assuming that the charmonium is absorbed with cross 

"The results for the LI sample were presented inE (obtained with a somewhat different analysis procedure). 
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section <to seeing in average < pL > amount of matter from its production point to its exit 
from the nucleus: a p A = (? Ae~ aabs<pL> (< pL >= (A — l)/2j ds T\{s) is obtained from the 
expansion of the Glauber formula). 

3) The widely used although not theory- motivated parametrization u p a = o Q A a . 



3 Results 



Table |l] shows the J ftp and ip' cross sections in the whole rapidity range (first column) and in 
four separate bins (columns 2-5). Due to the lack of space we show the cross sections obtained 
from the HI data sample only. Although the cross sections extracted from the LI data sample 
are systematically higher, it is the same ~ 5% difference both for the J/ip and the ip 1 ', which 
justifies the treatment of this discrepancy as a constant normalization factor. One should notice 
that since the systematic errors reflect the uncertainty in the absolute normalizations, they do 
not affect the shapes of the differential distributions in rapidity. 

Table 1: J/ip and cross sections per nucleon in nb (not corrected for the fifi branching ratio) for the integrated 
data and in four Ycm bins. The statistical and systematic errors are shown separately. 





-0.50< Y <0.50 


-0.50< Y <-0.25 


-0.25< Y < 0.0 


0.0< Y < 0.25 


0.25< Y < 0.5 




pBe 

pAl 

pCu 

pAg 

pW 


5.130±0.010±0.177 
4.868±0.008±0.228 
4.712±0.006±0.181 
4.403±0.005±0.148 
4.005±0.006±0.147 


1.202±0.007±0.042 
1.117±0.005±0.053 
1.173±0.004±0.045 
1.077±0.004±0.036 
0.945±0.004±0.035 


1.373±0.004±0.047 
1.304±0.003±0.061 
1.275±0.003±0.049 
1.190±0.002±0.040 
1.068±0.003±0.039 


1.352±0.004±0.047 
1.281±0.003±0.060 
1.209±0.003±0.047 
1.134±0.002±0.038 
1.047±0.003±0.038 


1.210±0.007±0.042 
1.145±0.005±0.054 
1.069±0.004±0.041 
1.016±0.004±0.034 
0.945±0.004±0.035 




pBe 

pAl 

pCu 

pAg 

pW 


.0886±.0021±.0032 
.0841±.0015±.0044 
.0773±.0011±.0032 
.0690±.0010±.0025 
.0611±.0010±.0024 


.0207±.0014±.0008 
.0194±.0010±.0009 
.0189±.0007±.0008 
.0151±.0006±.0005 
.0130±.0006±.0005 


.0231±.0010±.0008 
.0215±.0006±.0011 
.0210±.0005±.0008 
,0189±.0004±.0007 
.0161±.0004±.0006 


.0228±.0009±.0008 
.0228±.0006±.0012 
.0197±.0004±.0009 
,0182±.0004±.0007 
.0162±.0004±.0007 


.0217±.0013±.0009 
.0209±.0009±.0010 
.0173±.0007±.0009 
,0161±.0006±.0005 
.0158±.0007±.0006 



Fig. shows the J/ ip and ip' absorption parameters obtained from the global fit of HI and 
LI data samples by the models mentioned in Sec.^. For these joint fits the cross sections from 
the HI and LI samples were rescaled by the factors 2/(1 + R) and 2R/(1 + R), respectively 
(R = 0.957 is the result of the fit of the ratio of HI to LI values by a constant line). The results 
obtained separately from HI and LI samples are similar to the shown results. The large symbols 
in the center correspond to the fit of the data integrated over the —0.5 < y cm < 0.5 range. 
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Figure 1: a parameters and absorption cross sections for the J/ip and i/j' as a function of the Ycm obtained from 
the absolute cross sections fit: (Left) Glauber model; (Center) e~ <pL><Tabs ; (Right) A a parametrizations. The 
large symbols in the center correspond to fitting the data integrated over —0.5 < y cm < 0.5 . 
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Finally, Fig. [2] shows the difference between the a parameters and a a b s of the J/ip and 
ip' states extracted from fitting directly the ip' to J/ip cross section ratio (which cancels the 
uncertainties in the normalization). The full Glauber model is not applicable to cross sections 
ratios, so we show only the results of the e~ <pL>(7abs parametrization. 



Figure 2: Difference in the a parameters (Left) and in the absorption cross sections from the e~ <pL><Tabs 
parametrization for the J/ip and ?// states as a function of Yqm- Obtained from the direct fit of the ip' to 

J/ip cross sections ratio. 



4 Summary 

The production cross sections of J/ip and ip' in 450 GeV pA interactions was obtained and 
analyzed in the framework of charmonia nuclear absorption models. The J/ip nuclear absorption 
in the —0.5 < y cm < 0.5 range is consistent with previous NA50 results and yields a J '^ = 
0.927 ± 0.012 (for <j p a = cr A a parametrization). The fit with the full Glauber model yields 
a abt = 4-7 i 0.8 mb (or 4.3 ± 0.7 mb for the e~ <pL>Uabs parametrization). The corresponding 
fits for the ip' show a stronger absorption than for the J/ip, thus confirming the results of the 
E866 collaboration: of = 0.886 ±0.014 (with a^ bs = 7.7 ±0.8 and 6.7 ±0.8m6 for the Glauber 
and the e~ <pL>Uabs models respectively). The more precise fit of the ratio of ip 1 to J/ip cross 
sections yields a J ^ - of = 0.045 ± 0.007 and o^ hs - a J Jf = 2.6 ± 0.4 nb (for the e - <pL>a ^ 
parametrization). No statistically significant dependence of the absorption on the rapidity is 
seen, although there is an indication of a stronger absorption of the ip' at small rapidities. 
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